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ABSTRACT 



O ■ We investigate the dependence of the strength of galaxy clustering on intrinsic luminos- 

ity using the Anglo- Australian two degree field galaxy redshift survey (2dFGRS). The 
2dFGRS is over an order of magnitude larger than previous redshift surveys used to 
address this issue. We measure the projected two-point correlation function of galax- 
ies in a series of volume-limited samples. The projected correlation function is free 
from any distortion of the clustering pattern induced by peculiar motions and is well 
described by a power-law in pair separation over the range 0.1 < (r/h^ 1 Mpc) < 10. 
The clustering of L* (M^ — 51og 10 h = —19.7) galaxies in real space is well fit by a 
correlation length r$ = 4.9 ± 0.3 /i -1 Mpc and power-law slope 7 = 1.71 ± 0.06. The 
clustering amplitude increases slowly with absolute magnitude for galaxies fainter than 
M*, but rises more strongly at higher luminosities. At low luminosities, our results 
agree with measurements from the SSRS2 by Benoist et al. However, we find a weaker 
dependence of clustering strength on luminosity at the highest luminosities. The cor- 
relation function amplitude increases by a factor of 4.0 between M^, — 5 log 10 h = —18 
and —22.5, and the most luminous galaxies are 3.0 times more strongly clustered than 
L* galaxies. The power-law slope of the correlation function shows remarkably little 
variation for samples spanning a factor of 20 in luminosity. Our measurements are in 
very good agreement with the predictions of the hierarchical galaxy formation models 
of Benson et al. 

Key words: methods: statistical - methods: numerical - large-scale structure of 
Universe - galaxies: formation 



1 INTRODUCTION 

A major obstacle to be overcome by any successful theory of 
the formation of large scale structure is the problem of how 
galaxies trace the distribution of matter in the Universe. 



Measurements of differential galaxy clustering as a function 
of colour (Willmer et al. 1998), morphological type (Davis 
& Geller 1976; Iovino et al. 1993) and selection passband 
(Peacock 1997; Hoyle et al. 1999) imply the existence of 
biases between the distribution of galaxies and of mass. 
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A generic prediction of hierarchical structure formation 
models is that rarer objects should be more strongly clus- 
tered than average (Davis et al. 1985; White et al. 1987). 
Correspondingly, if more luminous galaxies are associated 
with more massive haloes, then these galaxies are expected 
to exhibit stronger clustering than the galaxy population as 
a whole (for the special case of bright galaxies at high red- 
shift, see for example Baugh et al. 1998; Governato et al. 
1998). However, the form of the dependence of the ampli- 
tude of galaxy clustering on luminosity remains controver- 
sial even after more than twenty years of constructing and 
analysing redshift surveys of the local Universe. In the liter- 
ature, claims of a dependence of galaxy clustering on lumi- 
nosity (e.g. Davis et al. 1988; Hamilton 1988; Maurogordato 
& Lachieze-Rey 1991; Park et al. 1994; Benoist et al. 1996; 
Willmer et al. 1998; Guzzo et al. 2000) have been made with 
similar regularity to claims of non-detections (e.g. Phillips 
& Shanks 1987; Hasegawa & Umemura 1993; Loveday et al. 
1995; Szapudi et al. 2000; Hawkins et al. 2001). Part of the 
reason for this disagreement is a mismatch in the range of 
luminosities and clustering length scales considered in ear- 
lier studies. However, the main problem with earlier work is 
the small size of the redshift surveys analysed, both in terms 
of volume and number of galaxies. With previous surveys, 
the dynamic range in luminosity for which clustering can 
be measured reliably is limited, particularly when volume- 
limited samples are used. Due to the small volumes probed, 
it has generally not been possible to compare the clustering 
of galaxies of different luminosity measured within the same 
volume. These results have generally been affected by sam- 
pling fluctuations that are difficult to quantify. This prob- 
lem is compounded by underestimation of the errors on the 
measured correlation functions and on the power-law fits 
traditionally employed in this subject. 

In this paper, we use the largest extant local survey, 
the Anglo- Australian two degree field galaxy redshift survey 
(hereafter 2dFGRS), to address the issue of how clustering 
depends upon galaxy luminosity. We describe the 2dFGRS 
and the construction of volume- limited samples in Section ^[ 
and our estimation of the correlation function is described in 
Section ^[ Our results for the real space correlation function 
are given in Section We compare our results with those 
from previous studies and with the predictions of simulations 
of hierarchical galaxy formation in Section 0. 



2 THE DATA 

2.1 The 2dFGRS sample 

The 2dFGRS is selected in the photometric b.j band from 
the APM galaxy survey (Maddox et al. 1990a,b; 1996) and 
its subsequent extensions (Maddox et al. in preparation). 
The survey is divided into two regions and covers approx- 
imately 2000 square degrees. The bulk of the solid angle 
of the survey is made up of two broad strips, one in the 
South Galactic Pole region (SGP) covering approximatively 
-37°.5 < S < -22°.5, 21 h 40 m < a < 3 h 30 m and the other in 
the direction of the North Galactic Pole (NGP), spanning 
-7°.5 < 8 < 2°5, 9 h 50 m < a < 14 h 50 m . In addition to these 
contiguous regions, there are a number of circular 2-degree 
fields scattered pseudo-randomly over the full extent of the 



low extinction regions of the southern APM galaxy survey. 
In this paper, we use the redshifts obtained prior to January 
2001, over 160 000 in total. As we are mainly interested in 
measuring clustering out to separations of order 20 Mpc, 
we do not include galaxies that lie in the random fields in 
our analysis. 

In order to select an optimal sample for the measure- 
ment of the two point correlation function, we apply a 
weighting scheme to objects in the 2dFGRS. A weight is 
assigned to each measured redshift based upon the redshift 
completeness mask, the construction of which is explained 
in Colless et al. (2001; see also Norberg et al. 2001, in prepa- 
ration). We require a relatively high completeness in a given 
direction on the sky, so that, in practice, our results are fairly 
insensitive to the precise details of the weighting scheme. Ex- 
cluding areas below our completeness threshold (which arise 
mainly as a result of the tiling strategy adopted to make op- 
timal use of telescope time, coupled with the fact that the 
survey is not yet finished), we estimate the effective solid 
angle used in the SGP region is ~ 420 n°, and in the NGP 

~ 190 n°. 

2.2 Constructing a volume-limited sample 

In this paper, we analyse a series of volume-limited subsam- 
ples drawn from the 2dFGRS. The advantage of this ap- 
proach is that the radial selection function is uniform, and 
the only variations in the space density of galaxies within 
each volume are due to clustering. By contrast, in a flux- 
limited survey, the galaxy number density is a strong func- 
tion of radial distance and this needs to be corrected for 
when measuring the clustering. The disadvantage of using 
a volume-limited sample is that a large number of galaxies 
in the flux-limited survey do not satisfy the selection cuts 
(which are explained below). This is a serious problem for 
previous surveys, but not for a survey the size of the 2dF- 
GRS. As we demonstrate in Section q, the volume-limited 
samples we analyse give robust clustering measurements and 
contain over an order of magnitude more galaxies than sim- 
ilar samples constructed from previous surveys (see Table 
!)• 

The construction of a volume-limited sample drawn 
from a flux-limited redshift survey requires a range of abso- 
lute magnitudes to be specified. Since a flux-limited survey 
has both bright and faint apparent magnitude limits, the 
selected range of absolute magnitudes requires that both 
a minimum (z m i n ) and a maximum (z max ) redshift cut be 
applied to the volume-limited sample. Thus, in principle, 
a galaxy included in the volume-limited sample could be 
displaced to any redshift between z m i n and z max and still re- 
main within the bright and faint apparent magnitude limits 
of the flux-limited survey. 

In order to estimate the absolute magnitude of 2dFGRS 
galaxies at redshift zero, it is necessary to apply corrections 
for band shifting (k— correction) and evolution in the stellar 
populations (e— correction). We adopt a global k + e correc- 
tion given by k + e = 0.03z/ (0.01 + z 4 ), which is a good 
fit to the correction calculated for the b.j selected ESO Slice 
Project survey using population synthesis models (see Fig. 
1 of Zucca et al. 1997). This form for the k + e correction 
gives consistent luminosity functions for the 2dFGRS when 
the survey is divided into redshift bins, indicating that it 
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Table 1. Properties of the combined NGP & SGP volume-limited sub-samples analysed. The second column gives the median magnitude 



of each sample. Columns 6 and 7 list the best fitting correlation length, ro, and power-law slope 7 of 
space, fitted over the range 0.5 < a/(h~ l Mpc) < 10. Column 8 gives the value of A(-y), defined by eqn 
value of 7. 
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adequately accounts for the degree of evolution in galaxy 
luminosity over the lookback time spanned by the survey 
(Norberg et al. 2001, in preparation). Our results are un- 
changed if we use the mean of the k corrections for different 
spectral types given by Madgwick et al. (2001, in prepara- 
tion). The values of z m i n and z max that define a volume- 
limited sample drawn from the 2dFGRS vary slightly with 
position on the sky. This is due to revisions made to the map 
of galactic extinction (Schlegel, Finkbeiner & Davis 1998) 
and to the CCD calibration of APM plate zero points since 
the definition of the original input catalogue. Throughout 
the paper, we adopt an Qo = 0.3, Ao = 0.7 cosmology to 
convert redshift into comoving distance. 



3 ESTIMATING THE TWO-POINT 
CORRELATION FUNCTION 

The galaxy correlation function is estimated on a two di- 
mensional grid of pair separations parallel (n) and perpen- 
dicular (a) to the line-of-sight. To estimate the mean density 
of pairs, a catalogue of unclustered points is generated with 
the same angular selection and (a m in, -Zmax) values as the 
data. The correlation function is estimated by 



DD - 2DR + RR 
RR ' 



(1) 



where DD, DR and RR are the suitably normalised number 
of weighted data-data, data-random and random-random 
pairs respectively in each bin (Landy & Szalay 1993). 

Contours of constant clustering amplitude in the red- 
shift space correlation function, f(<7, n), are distorted as a 
result of the peculiar motions of galaxies, as demonstrated 
for the 2dFGRS by Peacock et al. (2001). On small scales, 
random motions inside virialised structures elongate the 
constant-^ contours in the n direction, whereas on large 
scales, coherent flows flatten the contours. The latter ef- 
fect was measured clearly for the first time for galaxies us- 
ing the 2dFGRS (Peacock et al. 2001). The dependence of 
the redshift space correlation function on galaxy luminos- 
ity is analysed in a separate paper (Hawkins et al. 2001, in 



preparation). In this paper, to simplify the interpretation, 
we consider only clustering in real space, which we infer by 
projecting the measured correlation function along the line- 
of-sight. We compute a dimensionless quantity, H(<r) /a by 
integrating over the measured £(cr, n) grid (note that H(cr) 
is sometimes referred to as w(r p ) in the literature): 



'A°2 
a 



£(<r, 7t)d7r. 



(2) 



In practice, the integral converges by a pair separation of 
ty — 75 h~ 1 Mpc. The projected correlation function can, in 
turn, be written as an integral over the spherically averaged 
real space correlation function, £(r), 



a 



rdr 



(H 



(3) 



(Davis & Peebles 1983). If the real space correlation func- 
tion is a power-law (which is a reasonable approximation for 
APM galaxies out to separations around r ~ 10 h -1 Mpc, see 
e.g. Baugh 1996), then 



H(<x) _ fr y r(l/2)r([ 7 -l]/2) 

r(7/2) 



(4) 



where £(r) = (ro/r) 7 and ro is the correlation length. 

Previous studies have estimated the error on the mea- 
sured correlation function from the Poisson statistics of the 
pair counts in each bin (Peebles 1980) or by bootstrap re- 
sampling of the data (e.g. Benoist et al. 1996). Since we 
study a range of samples corresponding to different luminos- 
ity bins and also compare samples from different volumes, 
it is important to include an estimate of the sampling fluc- 
tuations in the error budget for the correlation function. 
This we derive from analysis of 22 mock 2dFGRS cata- 
logues constructed from the ACDM Hubble Volume dark 
matter simulation, in the manner explained by Baugh et al. 
(2001, in preparation; see also Cole et al. 1998). In order 
to mimic the clustering of the 2dFGRS, a biasing scheme is 
employed to select particles in the simulations with a prob- 
ability which is a function of the final dark matter density 
field, smoothed with a Gaussian filter (model 2 of Cole et al. 



© 0000 RAS, MNRAS 000, 000-000 



4 Norberg et al. 




0.5 1 

log 10 (a/h^Mpc) 

Figure 1. (a) The projected correlation function measured for 
galaxies in three different absolute magnitude bins in the same 
volume. The faintest sample contains 16 134 galaxies, the middle 
sample contains 6 186 galaxies and the brightest sample contains 
985 galaxies. For clarity, error bars are plotted only on the correla- 
tion function of galaxies with —18.5 > Mj, 7 — 51og 10 h > —19.5. 
(b) A comparison of the correlation function of galaxies in the 
same absolute magnitude bins but measured in different (al- 
though not completely independent) volumes. The heavy lines 
show results for galaxies with —18.5 > Mfcj — 51og 10 h > —19.5 
and the light lines show results for a brighter bin with —20.5 > 
Mb, — 51og 10 /i > —21.5. In each case, the dashed line shows 
the estimate from the optimal sample (see text) for the selected 
magnitude bin, whilst the solid line shows an estimate of the 
correlation function from the volume analysed in Fig. |l](a). For 
the —20.5 > Mf,j — 51og 1(1 /i > —21.5 magnitude bin, the opti- 
mal estimate is measured using 10 962 galaxies, which should be 
contrasted with the 985 galaxies used to make the measurement 
shown by the light solid line, in a volume defined by a broader 
magnitude bin. 



1998). The mock catalogues have the same clustering ampli- 
tude as galaxies in the flux-limited 2dFGRS and the same 
selection criteria that are applied to the data are used in 
the construction of the mock surveys. The clustering ampli- 
tude in the mocks is independent of luminosity. The error 
bars that we plot on correlation functions measured from 
the 2dFGRS are the rms found by averaging over the 22 
mock catalogues. 




0.2 0.4 0.6 0.8 1 

log^cj/tr'Mpc) 

Figure 2. The ratio of the projected correlation function of 
galaxies in different magnitude slices to the projected correlation 
function of galaxies with —19 > M bj — 51og 10 h > —20. Note that 
the ratio is plotted on a linear scale, whilst the pair separation 
is on a log scale. The solid line shows the ratio for galaxies with 
absolute magnitudes in the range —18 > Mf, — 51og 10 h > —19, 
the dotted line for —20 > M bj — 51og 10 h > —21 and the dashed 
line for —21 > Afj, — 51og 10 h > —22. For clarity, error bars have 
been omitted from the dotted line but these are comparable in 
size with those plotted on the solid line. 



4 RESULTS 

We first demonstrate the robustness of the approach of mea- 
suring the correlation function in volume-limited samples. 
Unless stated otherwise, we have added the pair counts in 
the NGP and SGP regions to compute correlation functions. 
In Fig. ^(a), we show the correlation function of galaxies 
in three disjoint absolute magnitude bins measured in the 
same volume. The sampling fluctuations are therefore virtu- 
ally the same for each sub-sample, although the number of 
galaxies varies between them. There is a clear difference in 
the clustering amplitude of galaxies in the brightest absolute 
magnitude bin. Next, we demonstrate that sampling fluctu- 
ations are not important in a survey the size of the 2dF- 
GRS. We show, in Fig. [l](b), the correlation function in two 
fixed absolute magnitude bins measured in different volume- 
limited sub-samples. Specifically, the dashed lines show the 
correlation function for the optimal volume-limited sample, 
appropriate to the selected absolute magnitude bin. Such a 
sample contains the maximum number of galaxies in that 
magnitude bin. The different estimates of the correlation 
function agree within the errors. 

We now focus attention on the series of volume-limited 
subsamples covering the range — 18 > M( JJ — 51og 10 /i > 
—22.5, whose characteristics are listed in Table 1. The 
shape and amplitude of the projected correlation function 
in a selection of these samples is compared in Fig. ^ with 
the correlation function of galaxies in the magnitude range 
— 19 > Mb r — 51og 10 /i > —20. The shape of the correla- 
tion function varies relatively little with the absolute mag- 
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Figure 3. (a) The correlation length in real space as a function of absolute magnitude. Results are shown for the SGP and NGP 
regions separately. The NGP points are plotted with an offset of 0.05 mag for clarity. Horizontal error bars on the SGP points indicate 
the absolute magnitude range of each bin, and each point is plotted at the bin centre. In both cases, the brightest data points are for 
galaxies in one magnitude wide bins. The solid line shows the predictions of the semi-analytic model of Benson et al. (2001), computed 
in a series of overlapping bins, each 0.5 magnitudes wide. The dotted curves show an estimate of the errors on this prediction, including 
the sample variance expected for a volume equal to that of the N-body simulation used, (b) The real space correlation length estimated 
combining pairs counts in the NGP and SGP (filled circles). The open symbols show a selection of recent data from other studies. The 
data for surveys selected in the B-band have been corrected to the bj band using the approximate relation M bj = Mg — 0.2. In order 
to compare samples defined by cumulative and differential magnitude bins, the data points are plotted at the median magnitude of each 
sample. 



nitude that defines the sample in contrast to the amplitude 
of the correlation function, which changes significantly for 
the brightest magnitude slice. Another view of this trend 
is given in Fig. 0(a) where we plot the real space correla- 
tion length as a function of absolute magnitude. The best 
fitting values of the correlation length, ro, and power-law 
slope 7, are determined by applying eqn. ^ to the mea- 
sured correlation function over the pair separation range 
0.5 < a/( h^ 1 Mpc) < 10 and carrying out a \ 2 minimisa- 
tion. This simple \ 2 approach will not, however, give reliable 
estimates of the errors on the fitted parameters due to the 
correlation between the estimates at differing pair separa- 
tions. We use the mock 2dFGRS catalogues to estimate the 
errors on the fitted parameters. In brief, the best fitting val- 



ues of ro and 7 are found for each mock individually, using 
the simple x 2 analysis. The estimated error bar is the rms 
scatter in the fitted parameters over the ensemble of mock 
catalogues. 

In Fig. [](a), we plot the correlation lengths for the NGP 
and SGP regions separately. These independent estimates 
are in excellent agreement with one another. The slope of the 
best fitting power-law correlation function, given in Table 
1, is similar for all the volume-limited samples considered. 
The clustering amplitude increases slowly with luminosity 
for galaxies fainter than M* (where M* = Mbj — 51og 10 h = 
— 19.7, as found by Folkes et al. 1999), but rises strongly 
at higher luminosities. The correlation function amplitude 
increases by a factor of 4.0 between M bj — 51og 10 h = — 18 
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0.5 



L/L* 

Figure 4. The variation of the relative bias as a function of lu- 
minosity, using the clustering of L* galaxies as a reference point 
(see text for definition). The 2dFGRS points are plotted at the 
median magnitude of each sample and the horizontal bars show 
the quartile magnitude range. The Benoist et al. (1996) points 
are taken from their Fig. 5 and are plotted at the median value of 
L/L* for each sample. Note that the error bars on the Benoist et 
al. points are obtained by averaging over correlated bins in pair 
separation. The curves shows parametric fits: the Benoist et al. 
measurements are well fitted by b/b* = 0.7 + 0.3L/L* (dashed 
line), whereas the 2dFGRS results suggest a more modest depen- 
dence on luminosity: b/b* = 0.85 + 0.15L/L* (solid line). 



and —22.5, and the most luminous galaxies are 3.0 times 
more strongly clustered than M* galaxies. 



5 DISCUSSION 

The volume-limited samples analysed in this paper contain 
over an order of magnitude more galaxies than previous 
studies of the dependence of clustering on galaxy luminos- 
ity, allowing a more accurate measurement of this effect than 
was possible before. The sheer volume covered by our sam- 
ples, 10 6 — 2 x 10 7 h~ 3 Mpc 3 , ensures that sampling fluctu- 
ations have little impact upon our results. 

We compare the 2dFGRS results with a selection of re- 
cent measurements taken from the literature since 1995 in 
Fig. [](b). To compare samples defined by cumulative and 
differential magnitude bins, we plot the datapoints at the 
median magnitude for the sample, as computed using the 
Schechter function parameters for the 2dFGRS (Folkes et al. 
1999). The horizontal bars plotted on selected points show 
the quartile range of the magnitude distribution in the sam- 
ple. Benoist et al. (1996) analysed quasi volume-limited sam- 
ples in the SGP region of the Southern Sky Redshift Survey 
2 (SSRS2), and found a sharp increase in the correlation 
length for galaxies brighter than Mb — 5 log 10 h — —20.5. 
The Benoist et al. correlation lengths are measured in red- 
shift space, although the authors report that a similar trend 
with luminosity is seen in real space. Willmer et al. (1998) 



re-analysed the SSRS2 South using different volume limits 
and also measured clustering in the SSRS2 North, present- 
ing fits for the correlation length in real and redshift space. 
Intriguingly, Willmer et al. find a larger correlation length 
in real space for galaxies with Mb — 51og 10 h ~ — 20 than 
Benoist et al. find in redshift space. Moreover, the clear 
disagreement between the results for the brightest galax- 
ies analysed in SSRS2 North and South suggests that sam- 
pling fluctuations are significant in a survey of this size and 
that the errors on these points have been underestimated 
(as demonstrated in Fig. 4 of Benson et al. 2001). Loveday 
et al. (1995) measured the clustering in real space by cross- 
correlating galaxies in the sparsely-sampled Stromlo/APM 
redshift survey with galaxies in the parent catalogue. Galax- 
ies were considered in three absolute magnitude bins. No 
difference was found between the clustering amplitude of 
L* and super-L* galaxies. However, the median magnitude 
for the most luminous sample considered by these authors 
is only 0.5 magnitudes brighter than M* . The increase in 
clustering amplitude with luminosity is connected with a 
change in the mix of morphological types with increasing 
luminosity. The mix of spectral types at the brightest abso- 
lute magnitudes is dominated by spectra characteristic of el- 
liptical galaxies, whereas spiral galaxies are more numerous 
around L* (Folkes et al. 1999; Cole et al. 2001; Madgwick et 
al. in preparation). The clustering of galaxies as a function 
of spectral type will be analysed in a separate paper. 

Our clustering results can be characterised in a concise 
way in terms of a relative bias parameter, b/b* , that gives the 
amplitude of the correlation function relative to that of L* 
galaxies (where M* = Mb, — 51og 10 h — —19.7). The relative 
bias between the correlation functions of galaxies of different 
luminosity is assumed to be constant for pair separations 
spanned by the ro values listed in Table 1 (see also Fig. 
2). The relative bias is then defined by b/b* — (Vo/?'o) 7/ ' 2 > 
where we take = 4.9 ± 0.3 ft" 1 Mpc from Table 1 and 
use 7 = 1.7. The 2dFGRS results are shown by the filled 
symbols in Fig. ^ and are well fitted by the relation b/b* — 
0.85 + 0.15L/L*. The 2dFGRS data suggest a significantly 
weaker dependence of the relative bias on luminosity than 
the Benoist et al. data, which follow the relation b/b* — 
0.7 + 0.3L/L* (Peacock et al. 2001). (The parametric fit 
to the Benoist et al. measurements was used by Peacock 
et al. 2001 to estimate the parameter /3 = Q 0M /b for L* 
galaxies in the 2dFGRS. Using the above fit to the 2dFGRS 
measurements changes the inferred value for f3 by less than 
la to (3 = 0.49 ± 0.08.) 

Hierarchical models of galaxy formation predict that 
bright galaxies should be more strongly clustered than faint 
galaxies (e.g. White et al. 1987; Kauffmann, Nusser & Stein- 
metz 1997). This generic prediction arises because bright 
galaxies are expected to occupy more massive dark matter 
haloes and these haloes are more strongly clustered than the 
overall distribution of dark matter. The trend of clustering 
amplitude with luminosity measured for 2dFGRS galaxies is 
in very good agreement with the predictions of a simulation 
of hierarchical galaxy formation taken from Fig. 4 of Benson 
et al. (2001), reproduced as the solid line in Fig. ^(a). In the 
Benson et al. semi-analytic model, the input parameters are 
set in order to reproduce a subset of local galaxy data, with 
most emphasis given to the field galaxy luminosity function 
(see Cole et al. 2000). No reference is made to clustering 
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data in setting the model parameters. In a ACDM cosmol- 
ogy, Benson et al. (2000a,b) find excellent agreement with 
the real space correlation function measured for galaxies in 
the APM survey by Baugh (1996). It is remarkable that the 
same model, without any readjustment of parameters, also 
reproduces the dependence of clustering amplitude on lumi- 
nosity exhibited by the 2dFGRS in Fig. |(a). 
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